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In a mini workshop on B Physics which were held at Dalian, China, a number of lectures was given. 
Around 40 participants joined the workshop. In this workshop, people discussed the heavy quark 
effective field theory, light cone QCD sum rules and light cone wave functions. Much discussion 
was devoted to the currently popular QCD factorization and perturbative QCD approaches in the 
hadronic B decays. The physics of CP violation and physics beyond standard model were also 
discussed. 

I. OUTLINE 

The standard model (SM) has been well tested in the gauge sector, the unsolved and unclear problems in the SM 
are mainly concerned in the Yukawa sector which is strongly related to flavor physics, such as: origin and mechanism 
of CP violation, origin of the quark and lepton masses as well as their mixing. It involves thirteen parameters whose 
origin are all unknown. Therefore, precisely extracting those parameters and testing CP violation mechanism as well 
as probing new physics become hot topics in flavor physics. In fact, flavor physics has already indicated the existence 
of new physics. After forty years of discovery for indirect CP violation, direct CP violation in kaon decays has well 
been established. Theoretical predictions p| and experimental measurements 0, Q are now consistent each other. 
Exclusive semi-leptonic and inclusive B decays play a crucial role for extracting two important parameters Vd, and 
V u b in the CKM matrix elements. Rare B decays and direct CP violations are also of great importance in determining 
weak phase angles of the unitarity triangle and testing the Kobayashi-Maskawa (KM) mechanism |^| in SM as well 
as probing new physics. 

The currently running B-Factories at KEK and SLAC provide more and more experimental data on B physics. 
The future LHCb, BTeV and super-B factory will give much more data. More Precise Experimental Data and more 
precision theoretical prediction will encourage precision test of SM at B sector. 

In this workshop (next section), Y.L. Wu talks about how the heavy quark effective field theory provides a powerful 
tool for studying beauty physics, what is the implication of beauty physics. Some attentions are also paid to the more 
precise extraction of V c b and V u b and large direct CP violation in charmless B decays. 

In section lITn T. Huang gives a brief review on the light-cone wave function(LCWF) in QCD. The definition of the 
light-cone wave function in QCD and distribution amplitude and its asymptotic form are reviewed. General properties 
of the light-cone wave function and higher helicity components are discussed. C.F. Qiao discusses about the derivation 
of B meson wave function. X.H. Wu talks about the twist 3 wave function of pion in QCD sum rules. The 1/rrib 
power suppressed effects are discussed by Z.H. Li. 

Y.D. Yang introduces the QCD factorization approach and discusses the application in hadronic B decays in section 
IIVI K.T. Chao intensively studies the B decays to charmonium final states in QCD factorization in section^ X.Q. Yu 
and CD. Lu present the formalism and application of perturbative QCD approach, which is based on kx factorization 
in section IVTl Intensive discussions are induced for the comparison and underlining theory of these two approaches. 

The CKM angle determination is one of the main talks for the B factories. Z.J. Xiao discusses a number of channels 
which are useful for these angles a, j3 and 7 in section IVTT1 In section IVlIII C.H. Chang discuss the exact solution of 
Bethe-Salpeter (BS) equation. 

The new physics picture is always a hot topic. X.G. He and C.S. Huang give the SUSY contributions in B — > K*j 
decay and B — > (j)Ks and B — ► rj'Ks decays in section llXl These decays are more sensitive to new physics than some 
others. Z.J. Xiao makes a systematic study for the new physics contributions to the charmless two-body hadronic 
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decays of B and B s meson induced by the neutral- and charged-Higgs penguin diagrams in model III: the third type 
of the two-Higgs-doublet models and technicolor models. 

II. LCQCD/HQEFT FROM FULL QCD 

by Y.L. Wu 

Heavyquark effective field theory of QCD which was first explored in [j| and detailed developed recently in a serious 
papers [y, 0, M> M, OH LUl U4 U& G_J HE Un LLIl provides a premising and systematic tool in correctly evaluating the 
hadronic matrix elements of heavy quarks and extracting the CKM matrix elements V c b and V u b from B decays via 
heavy quark expansion (HQE). HQEFT is a theoretical framework derived directly from QCD but explicitly displays 
the heavy quark symmetry in the infinite mass limit ttlq — > oo and symmetry breaking corrections for finite mass 
case in the real world. In fact, HQEFT has been shown to be as a large component QCD 0, |^. At the leading 
order, it coincides with the usual heavy quark effective theory(HQET) jl8j which is constructed based on the heavy 
quark symmetry (HQS) 0] in the infinite mass limit [20|. The differences between HQEFT of QCD and HQET arise 
from the sub- leading terms in the l/rag expansion. This is because in the construction of HQET the particle and 
antiparticle components were separately treated based on the assumption that the particle number and antiparticle 
number are conserved in the effective Lagrangian plj . However, such an assumption is only valid in the infinite 
mass limit. Note that the particle number and antiparticle number are always conserved in the transition matrix, 
which is independent of heavy quark limit. Therefore the particle-antiparticle coupling terms which correspond to the 
pair creation and annihilation terms in full QCD were inappropriately dropped away in the HQET though they are 
suppressed by 1/toq but only vanish in infinite mass limit. In HQEFT of QCD 5, [15J, all contributions of the field 
components, large and small, 'particle' and 'antiparticle', have carefully been treated in the effective Lagrangian, so 
that the resulting effective Lagrangian form the basis for a "complete" effective field theory of heavy quarks. It is seen 
that the particle-antiparticle coupling terms are suppressed by 1/toq and they are decoupled in the infinite mass limits. 
Their physics effects at sub-leading order have been found to be significant in some cases [sL Fil IM 1^. Hoi ITU IT3L Hif . 
Therefore, to consider the finite quark mass correction precisely, it is necessary to include the contributions from the 
components of the anti-quark fields. Note that the Wilson coefficient functions describe the perturbative effects in full 
QCD, HQEFT of QCD treats the non-perturbation effects below the energy scale nib- So that the anti-quark effects 
in HQEFT of QCD cannot be attributed to the Wilson coefficient functions. As the anti-quark effects appear from 
the sub- leading order of 1/toq expansion, it should not surprise that at the leading order the resulting anomalous 
dimensions from both full QCD and HQEFT of QCD are all the same. 

A complete description for the theoretical framework of LCQCD/HQEFT is going to be presented in a longer paper 
with the following contents: 

1. Introduction 

2. Effective Lagrangian of LCQCD/HQEFT from QCD 

3. Lorentz Invariance of LCQCD/HQEFT 

4. Quantization of LCQCD/HQEFT 

4.1 Quantum Generators of Poincare Group 

4.2 Anti-commutations and Super selection Rule 

4.3 Hilbert and Fock Space of LCQCD/HQEFT 

4.4 Propagator in LCQCD/HQEFT 

4.5 Discrete Symmetries in LCQCD/HQEFT 

5. Basic Framework of LCQCD/HQEFT 

5.1 Feynman Rules in LCQCD/HQEFT 

5.2 Finite Mass Corrections 

5.3 Trivialization of Gluon Couplings and Decouple Theorem 

5.4 Renormalization in LCQCD/HQEFT and Wilson Loop 

5.5 Current Operators in LCQCD/HQEFT 

6. New Symmetries in LCQCD/HQEFT 

7. Demonstration for the Trace Formula of Transition Matrix Elements and Universal Isgur-Wise Function 

8. Simple applications 

Hereafter is a summary of Important Effects in LCQCD/HQEFT: 

• Automatically reparameterization invariance v — > v' with v 2 = v 2 = 1. 

• Automatically with vanishing 1/toq corrections at Zero- Recoil (both for B — > D*lvi & B — > Dlvi ) 
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• Largely reduce the numbers of form factors at higher Order of 1/ tjiq and more relations are obtained 

• Nicely keep scaling law of decay constants 



F -I 1 + O(sb)] (~ 10% for b . ~ 30% for c ) 

e as follov 

2dMff2)} 



The explicit forms and results were found to be as follows 
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F = 0.38 ±0.06 GeV 3/2 ; 
The numerical results are 



gi = 0.46 ± 0.12 GeV; g 2 = -0.06 ± 0.02 GeV 



f B (m h ) = 0.196 ± 0.044 GeV, 
f D (m c ) = 0.298 ± 0.109 GeV, 



f B , (m h ) = 0.206 ± 0.039 GeV, 
f D , (m c ) = 0.354 ± 0.090 GeV. 



• Reliably understand the lifetime differences of bottom hadrons B d , B s , A b . 



r(B°) 
r(B-) 



0.79 ±0.01 



1.08 ±0.05. 



0.96 ±0.06 



T(B d ) T{B d ) 

B S l{B°) = (10.51 ± 0.37), B SL {A b ) = (11.13 ± 0.26)% 
n c {B a ) = 1.20 ± 0.03, B T (B°) = (2.6 ± 0.1)% 



Consistently extraction of V c b and V u b 



0.0395 ± O.OOllexp ± 0.0019 th 
0.0434 ± 0.0041 exp ± 0.0020 th 
0.0394 ± 0.0010exp ± 0.0014 th 
0.0402 ± 0.0014 exp ± 0.0017 th 



V cb \ 

V cb | 

v cb \ 

V cb\ 

V ub \ LO = (3.4 ± 0.5 exp ± 0.5 th ) x lO" 3 
V ub \ LO = (3.7 ± 0.6 exp ± 0.7 th ) x 10~ 3 
V ub \ NLO = (3.2 ± 0.5 exp ± 0.2 th ) x 10" 3 
V ub \ = (3.48 ± 0.62 exp ± 0.11 th ) x lO" 3 



from B 

from B 

from B 

{Average) 
from B - 
from B - 
from B 
from B - 



D*lv, 
Dlu, 
X r lv, 



0(l/m 2 Q ) 
0(l/m|) 
0(l/ml) 



ply, 
-> ttIv, 

>XJv 0{l/m%) 



III. LIGHT-CONE SUM RULES AND LIGHT-CONE WAVE FUNCTION IN QCD 



by T. Huang, Z.H. Li, C.F. Qiao and X.H. Wu 

A brief review on the light-cone wave function in QCD is presented.lt consists of the following nine sections: 1) the 
definition of the light-cone wave function in QCD; 2) distribution amplitude(DA) and its asymptotic form; 3) general 
properties of the light-cone wave function; 4) Melosh transformation and higher helicity components; 5) DA moments 
in the QCD sum rules; 6) model approach to LCWF; 7) Conformal expansion; 8) Twist-3 distribution amplitude of 
the pion; 9) perspective on the light-cone wave function. 

As well known, the light cone formalism provides a convenient framework for the relativistic description of hadrons 
in terms of quark and gluon degrees of freedom, and the application of perturbative QCD (PQCD) to exclusive 
processes has mainly been developed in this formalism [24j].In the PQCD theory, the hadronic distribution amplitudes 
and structures which enter exclusive and inclusive processes via the factorization theorems at high momentum transfer 
can be determined by the hadronic wave functions, and they are the underlying links between hadronic phenomena 
in QCD at large distances (nonperturbative) and small distance (perturbative). A light-cone (LC) wave function is 
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defined at the LC time r = x + = x° + x 3 in physical light cone gauge A + = A + A 3 = 0, which is conjugate to the 
LC Hamiltonian Hlc = P = P° — P 3 [25J.A stationary solution |\P(t)) has a LC Fock state expansion and the LC 
wave functions, H! n (xi , kj_i , Ai), are the amplitudes to find n partons (quarks, anti-quarks and gluons) with momenta 
ki in a pion of momentum P. Hence it is important to develop methods which specify the light-cone wave functions 
of the hadron. 

The distribution amplitude &(x,Q) for the leading twist satisfies a QCD evolution equation [24j.As Q 2 goes to 
infinity, $ (x , Q) is dominated mainly by the first term in the expansion, i.e for the twist-2 DA of the meson, we have 
— > a^Jx(l — x). However the solution of the evolution equation can be determined completely only as if the 
initial condition of Qo) is given. 

From the relation between the wave functions and measurable quantities we can get the general properties of the 
hadronic wave functions |25|. For example, in the pion case two important constraints on the valence state wave 
function ^ q q(x, kj_) (for the Ai + A2 = components) have been derived from it — * fiu and ir° — > 77 decay amplitude. 
The quark transverse momentum of the valence state in the pion (kj_) g q should be larger than (k^)^. P q q is the 
probability of finding the qq Fock state in the pion, and it should be less than unity, the charged pion radius for the 
lowest valence quark state (r^ + ) 9 ' should be lower than the experimental value for (r 2 + ), where the latter one should 
contain all the contributions of the higher Fock stats. 

In order to obtain the LC spin wave function of the pion, we transform the ordinary instant-form SU(6) quark 
model wave function into a LC one. The instant-form spin states \J,s)t and the LC-form spin states \J,X)p are 
related by | J, \}p = ^2 S U^ X \J, s)t ■ This rotation is called the Melosh rotation for spin- 1/2 particles. Applying the 
Melosh rotation we can obtain the spin wave function of the pion in the infinite-momentum frame |27| . There are 
two higher hclicity components (Ai + A2 = ±1) in the expression of the LC spin wave function besides the ordinary 
helicity components (Ai + A2 = 0). These higher helicity components will make contributions to the exclusive processes 
although they have the power suppression behavior. 

Applying the QCD sum rules to the distribution amplitude,we can get the first three moments of the initial 
DA, $(2;, Qo) [2^|. From these moments we can not determine the initial &(x,Qo) and only know that the initial 
&{x, Qo) is different from the asymptotic form. 

For the space wave function, one can take the Brodsky-Huang-Lepage (BHL) prescription which is obtained by 
connecting the equal-time wave function in the rest frame and the wave function in the infinite momentum frame 

(25j. We obtain the LC space wave function of the pion, ifiBUL = A exp [— 8 ^-t^^_ - j ] [25j . This form provides a more 
reasonable Ut dependence in the space wave function and is consistent with the experimental data |28j . 

Conformal expansion allows one to separate transverse and longitudinal variables in the distribution amplitude. one 

1/2 3 /2 . P^l 

can expand DA in terms of Gegenbauer polynomials C n and C„ according to its conformal spin [2Jj . This symmetry 
helps us to make a model for the light-cone wave function of the hadron. 

Recently, we do not use the equation of motion for the quarks inside the pion and apply the QCD sum rules 
to calculate the pion twist-3 DA $p(^) and get its approximate expression and its expansion coefficients can be 
determined by its three moments approximately |3Cj . Based on the moment calculation we suggest a model for the 
twist-3 wave function with kx dependence |3l|. 

In Summary, at present although we can not solve light-cone wave function of the hadron completely due to its 
non-perturbative dynamics, one has got much progress by studying its general properties, symmetry, its moments from 
non-perturbative calculation and phenomenological model. These studies will contribute to a better understanding 
on the light-cone wave functions and approach their realistic solutions. 

A. Twist-3 Distribute Amplitude of the Pion in QCD Sum Rules 

We apply the background field method to calculate the moments of the pion two-particle twist-3 distribution 
amplitude (DA) </> p (£) in QCD sum rules. The QCD sum rules were used to study the leading-twist distribution 
amp litude of the pion at the first time in pfij . The pion twist-3 distribution amplitudes have been studied in Ref. 
|32l l33| in the chiral limit, based on the techniques of nonlocal product expansion and conformal expansion and 
including corrections in the meson-mass. However they employ the equations of motion of on-shell quarks in the 
meson to get two relations among two-particles twist-3 distribution amplitudes of the pion, (j)p(£.) an d 0«t(£)j an d 
three-particle twist-3 distribution amplitude 03„-(ai) of the pion. The question is whether one can use the equation of 
motion of the quark inside the meson since the quarks are not on-shell. In our paper [3(J, we do not apply the quark 
equation of motion and calculate the moments of the twist-3 distribution amplitudes <p p (£) directly in QCD sum rules 
approach. 
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For calculating the moments of 4>p{£) , we take as usual the two-point correlation functions: 

d*xe^(0 | T{ji 2n) (x,z)jP f (0)} | 0) = (z- q f n I^ Q (q 2 ), 
where the currents are defined as 

jf n) {x 1 z)=d{x) l5 {iz-^D) 2n u{x) 

withD = — i2gA®T a . To get the sum rules for moments, we employ dispersion relation for Ip 7l '°(Q 2 ) and isolate 
the pole term of the lowest pseudoscalar pion. With condensates up to dimension-6 and the perturbative contribution 
part to the lowest order, we obtain the sum rule for moments of </> p (£): 



M 4 



(£ 2n ) ■ = — — 



\/M 2 



8tt 2 2n 



+ 1 V M 2 J 



V j A/ 4 24 M 4 81 V V " M 6 

with to = (m„ + m<j)/2, and Af is a Borel parameter. For the condensate parameters, we take as usual. And the 
definition of the moments are given by 



2n \ 



The key point to avoid the concept of on-shell equations of motion is that the introduced mjj appears in the sum 
rule of moments. And they can be obtained through requiring normalization of the zeroth moments((£p) = 1). The 
advantage is that our results do not rely on the validness of the equation of motion and at the same time provide a 
parameter which can be compared with the corresponding one used in phcnomcnological analysis. With these points 
kept in mind, one can show that the results for the first three moments are(M 2 = 1.5 — 2 GeV 2 ): 



s w (GeV 2 ) 
1.7 
1.6 
1.5 



mg(GeV) 
1.24- 1.36 
1.19-1.30 
1.14-1.24 



0.340 - 0.356 
0.340 - 0.359 
0.341 - 0.361 



(£ 4 ) 

0.167-0.210, 
0.164-0.211, 
0.160-0.212. 



It can be seen that, the parameter is smaller than m^j{m u + md) — 1-78 — 3.92 GeV which was used in, e.g., j^] 
when the equation of motion of on-shell quarks are employed. We take the first three moments into consideration in 
our analysis. Expanding twist-3 distribution amplitudes 4> P {£,) °f the pion in Gegenbauer polynomials to first three 
terms, one can obtain its approximate form, 



^(0 = 1 + 0.137 C 2 1/2 (£) 



0.721 C 4 1/2 (£)- 



B. B-Meson Wave Functions 



The development of the heavy quark effective theory(HQET) has led to much progress in the theoretical under- 
standing of the properties of hadrons. With the progress in both theory and experiment, nowadays, B physics becomes 
one of the most active research areas in high energy physics. Many B meson exclusive decay processes turn out to 
be calculable systematically in the frameworks of newly developed factorization formalisms, the so-called PQCD ap- 
proach [34| and QCD Factorization approach j35|. In all of these calculations based on the factorization approaches, 
the light-cone distribution amplitudes of the participating mesons, which express the nonperturbative long-distance 
contributions in the factorized amplitudes, play an important role in making reliable predictions. 

By definition, the light-cone distribution amplitudes are given by the light-cone wave functions at zero transverse 
separation of the constituents. Following j3|| it is 



\q(z)Th v (0)\B(p)} 



if B M 



Tr 



(t) 



2t 



(1) 
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Eq. Q is the most general parametrization of the two-particle light-cone matrix element compatible with Lorentz 
invariance and the heavy-quark limit. The three-particle matrix element can be generally expressed as |37| : 



(0\q(z)gG^(uz)z v Th v (0)\B(p)) — — /bMTi 



75 r 



1+ 



K^-*7m) (*^(t,lt)-*v(*,u)) -ia„ u z u ^ v (t,u)- z„X A (t,u) + -± /Y A (t,u) 



(2) 



By virtue of the equations of motion, one can get a set of relations between distribution amplitudes. In the approxi- 
mation that the three-particle amplitudes are set to be zero, the Wandzura-Wilczek approximation, the system of the 
obtained differential equations is simplified, and the solutions for <p± can be obtained. In the momentum space, the 
analytic solution can be expressed explicitly as |37j 



ip+(u,k T ) 



OJ 

27rA 2 ' 



(w)0(2A - w)6 (k% - w(2A - to)) 



(3) 



Tjj-(uj,k T ) = ttt 6(u)d(2A - u)S - w(2A - u)) (A) 

The results (01 and J3J give exact description of the valence Fock components of the B meson wave functions in 
the heavy-quark limit, and represent their transverse momentum dependence explicitly. These results show that the 
dynamics within the two-particle Fock states is determined solely in terms of a single nonperturbative parameter A. 

It should be noted that the Wandzura-Wilczek approximation is not equivalent to the free field approximation. The 
leading Fock-states, which correspond to the twist-2 contribution in the case of light meson wave functions, carries 
the effect from the QCD interaction. 

There has been indication that, in the heavy-light quark systems, the higher Fock states could play important 
roles even in the leading twist level [3?} . This would suggest that the shape of the wave functions as function of 
momenta and their quantitative role in the phenomenological applications would be modified when including the 
higher Fock states. However, the direct applications of our wave functions to phenomenological processes indicate 
that the Wandzura-Wilczek approximation does not induce obvious conflicts with experimental data |38( 

In contrast to many model assumptions, which show strong damping at large \zt\ a $ ~ exp f— K 2 z^/2), our 
wave functions J3J and 10} have slow-damping with oscillatory behavior as ip±(u>, — z%) ~ cos(|zr| \A->(2A — u) — 
ir / A) / yf\zr\. On the other hand, heavy-quark symmetry guarantees that the solution in our work provides complete 
description of the light-cone valence Fock wave functions for the B* mesons and also for the D, D* mesons in the 
heavy-quark limit. Moreover, the investigations on beyond Wandzura-Wilczek limit are in progress. 



C. 1/rrib power-suppressed effects in two-body B decays 

A careful investigation of 1/ mj power suppressed effects is crucial for a good understanding of two-body hadronic 
decays of B mesons. It is instructive to systematically investigate the 1/mb suppressed effects from either soft or hard 
gluon exchanges, in addition to annihilation topology. Focusing on penguin-dominated B — > Kir and color-suppressed 
B a — > D a ir° decays, we estimate such power-suppressed corrections which appear possibly in the framework of QCD 
light-cone sum rules (LCSR) . 

In the former case, the soft exchanges could be induced by the emission topology or by the chromo-magnetic diploe 
operator Os g and penguin topology, while the hard gluon exchanges may occur through the Os g operator or penguin 
contractions. The behaviors of such contributions in the heavy quark limit nib - * oo comply with the following power 
counting: (1) The soft effects due to the emission topology are of (2) The penguin diagrams make no 

contributions up to l/r?i| order. (3) The Os g operator supplies a hard effect of 0(l/nib) and a soft effect of 0(1/ ml). 
Despite being formally suppressed by 0(l/ml) compared with the leading-order factorizable amplitudes, the soft 
correction from Og g is free of a s suppression and numerically comparable with the 0(a s ) hard part. (4) The part 
with quark condensate has a chiral enhancement factor r£ via the PCAC relation, suffering formally from 0(a s ) 
and 0(l/rrib) double suppression but numerically turning out to be a large effect. However, it has a counterpart in 
the QCD factorization and thus is not included in our calculation to avoid a double counting. For the weak phase 
7(= lmV* b ) ranging from 40° to 80°, we find that because of the annihilation and soft and hard exchanges a numerical 
increase of (20 — 30)% is expected in the branching ratios up to 0(l/ml) corrections; the resultant soft and hard 
corrections are less important than the annihilation contributions and amount only to a level of 10%. Possible sources 
of uncertainties are discussed in some details. 
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For the case of B° — > D°ir°, naive factorization assumption breaks down and non-factorizable contributions are 
leading. Then the 1 /mi, power-suppressed soft exchange between the emitted heavy- light quark pair and Btt system 
is expected to be much more important than in the former case. The resulting correction to the decay amplitude 
is found to be comparable numerically with the corresponding factorizable piece, estimated at about (50 — 110)% 
of the latter, and renders the relevant parameter 02 receive a positive number correction. This indicates that such 
power-suppressed effect would be crucial to our phenomenological understanding of the B° — ► D°ir° decay. 

IV. QCD IMPROVED FACTORIZATION APPROACH FOR B DECAYS 

by Y.D. Yang 

It is of great interest and importance to investigate the decays of B mesons to charmless final states to study the 
weak interactions and CP violation. In past years, we have witnessed many experimental and theoretical progresses 
in the study of B physics with the observations of many B charmless decay processes and improvements of the theory 
of B decays. 

Most of the theoretical studies of B decays to pseudoscalar and vector final states are based on the popular Naive 
Factorization approach |39|. As it was pointed out years ago in Ref . |4d| . the dominant contribution in B decays 
comes from the so-called Feynman mechanism, where the energetic quark created in the weak decay picks up the 
soft spectator softly and carries nearly all of the final-state meson's momentum. It is also shown that Pion form 
factor in QCD at intermediate energy scale is dominated by Feynman mechanism 0,E2>E3- iFrom this point, wc 
can understand why the naive factorization approach have worked well for B and D decays, and the many existing 
predictions for B decays based on naive factorization and spectator ansatz do have taken in the dominant physics 
effects although there are shortcomings. However, with the many new data available from CLEO and an abundance 
of data to arrive within few years from the B factories BaBar and Belle, it is demanded highly to go beyond the naive 
factorization approach. 

Recently, Beneke et.al., have formed an interesting QCD factorization formula for B exclusive non-leptonic de- 
cays |35l [y|- The factorization formula incorporates elements of the naive factorization approaches leading con- 
tribution) and the hard-scattering approaches sub-leading corrections), which allows us to calculate systematically 
radiative (sub-leading non-factorizable) corrections to naive factorization for B exclusive non-leptonic decays. An im- 
portant product of the formula is that the strong final-state interaction phases are calculable from the first principle 
which arise from the hard-scattering kernel and hence process dependent. The strong phases are very important for 
studying CP violation in B decays. Detail proofs and arguments could be found in [35j. 

The amplitude of B decays to two light mesons, say M\ and M2, is obtained through the hadronic matrix element 
(Mi(pi)M2(p2)\Oi\B(p)), here M\ denotes the final meson that picks up the light spectator quark in the B meson, 
and M2 is the another meson which is composed of the quarks produced from the weak decay point of b quark. Since 
the quark pair, forming M2, is ejected from the decay point of b quark carrying the large energy of order of m&, soft 
gluons with the momentum of order of Aqcd decouple from it at leading order of Aqq d / f fnb i n the heavy quark limit. 
As a consequence any interaction between the quarks of M2 and the quarks out of M2 is hard at leading power in the 
heavy quark expansion. On the other hand, the light spectator quark carries the momentum of the order of Aqcd, 
and is softly transferred into Mi unless it undergoes a hard interaction. Any soft interaction between the spectator 
quark and other constituents in B and M\ can be put into the transition form factor of B — > M\ . The non-factorizable 
contribution to B — ► M\M.2 can be calculated through the diagrams in Fig.l. 

It is worth to note that the shortcomings in the "generalized factorization" are resolved in this framework. Non- 
factorizable effects are calculated in a rigorous way here instead of being parameterized by effective color number. 
Since the hard scattering kernels are convoluted with the light cone DAs of the mesons, gluon virtuality k 2 = xm\ 
in the penguin diagram Fig.l.e has well defined meaning and leaves no ambiguity as to the value of fc 2 , which has 
usually been treated as a free phenomenological parameter in the estimations of the strong phase generated though 
the BSS mechanism |4^. So that CP asymmetries could be predicted soundly. 

V. B MESON EXCLUSIVE DECAYS TO CHARMONIUM IN QCD FACTORIZATION 

by K.T. Chao 

Exclusive i?-meson decays into charmonium are interesting, since these decays e.g. B — > J /ipK are regarded as 
the golden channels for CP violation studies, and these decays provide useful information towards the understanding 
of perturbative and nonperturbativc QCD. These decays are color-suppressed decays, and involve two heavy quark 



8 




M 




->-W, /0O5» 




B 


Mi 


M 




Mi 



(f) 





Figure 1: 



FIG. 1: Feynman diagrams contributed to non-factorizable contributions. 



energy scales, mj, and m c , therefore are more subtle in theoretical studies. While experimentally CLEO, BaBar and 
Belle Collaborations have provided many measurements on the B meson exclusive decays to charmonium [4^ | such as 
J/i>,ip' ,T]c,v'ciXcOi an d Xcii theoretical studies for those decays are still limited. 

In the QCD factorization approach [351 Eij. it is argued that because the size of charmonium is small (~l/a s mj,) 
and its overlap with the (B,K) system is negligible in the heavy charm quark limit, QCD factorization method can 
be used for B — > J/ipK and other charmonium states, where charmonium is described by the color-singlet cc pair. 
Indeed, explicit calculations for B — > J /ijjK within the QCD factorization approach 0] show that the nonfactorizable 
vertex contribution is infrared safe and the spectator contribution is perturbatively calculable at twist-2 order, but the 
theoretical branching ratio is much smaller than the experimental data. Aside from the B — > J/ijjK decay |47| , here 
we pre sent some additional results for B exclusive decays to other S-wave, P-wave, and D-wave charmonium states 
|48l l49t l50L l51| in QCD factorization. Problem of infrared divergences in B decays to P-wave charmonium (B ^XcoK : 
B — > Xc2K, and B — > h c K) and problem of S-D mixing in B decay to D-wave charmonium (B — > ^(3770)-ft')are 
emphasized. 



A. B -^rj c K decay 

Studies of B meson decays into pseudoscalar charmonium states B -^rj c K and B -^r] c (2S)K 4^ m QCD factor- 
ization show that the nonfactorizable corrections to naive factorization are infrared safe at leading-twist order. The 
spectator interactions arising from the kaon twist-3 effects are formally power-suppressed but chirally and logarith- 
mically enhanced. An important improvement by including the 0{a s ) corrections is the large cancellation of the 
renormalization scale /x dependence of the decay amplitude. However, the calculated decay rates are smaller than 
experimental data by a factor of 8-10. On the other hand, it is found that for B meson decays to J/tp,tp' ,rj c and 
r]' c , the predicted relative decay rates of these four states are approximately compatible with experimental data. For 
instance, the theoretical value 

Br(j?° -> Vc K°) _^_ 2 . FiKJ 2 ml- ml a (Jjh_) 2 ^ 83 - 1 2 (S) 



is compatible with the experimental value 

Br( J B° -> 7] C K°) 



Br(B° -► J/iPK°) j 



1.0. (6) 



, 2 



Here F\ is the B —t K form factor, and the ratio of the squared decay constants {jjj^Y ranges from 0.75 to 1.1 in 
various models. 
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B. B ~^XcoK, B —*xczK, B ^h c K decays 

Studies of B meson decays to P-wave charmonium states B — ► XcoK, B — > x C 2 AT, and B — > h c K 0, challenge the 
applicability of QCD factorization for charmonium because of the appearance of non-vanishing infrared divergences 
in these decays. 

For these decays, there are no leading order contributions from the V —A currents in the effective 4- qua rk lagrangian, 
because these P-wave charmonium states do not couple to the V — A currents. According to [3{J 0] in QCD 
factorization all nonfactorizable corrections are due to the vertex corrections and spectator corrections, and other 
corrections are factorized into the physical form factors and meson wave functions. Taking nonfactorizable corrections 
into account, the decay amplitude for B —>Xc.jK(J — 0, 2) can be written as 



iM = — = 

V2 



VcbV^d-VtbV^Ci + Ce) xA, (7) 
where Ci are the well known Wilson coefficients and A is given by 



*6Ki(0) * S C F( 4n 2 f K f B . , 

r-rr ■ XriTT [*l-fl + Kf in)- ( 8 ) 



iwN c V x J ' 1 N c 

Here R[ (0) is the derivative of the radial wave function at the origin of P-wave charmonia , N c is the number of colors, 
Cf = (AT? — 1)/(2N C ). The function // is calculated from vertex corrections, and /// is calculated from spectator 
corrections. For 3 Po charmonium state Xco, the infrared divergence from the vertex correction is found to be 

8mfcz(l — z + In z) A 2 

fi = -== ln( — ^) + finite terms, (9) 

(1 - z)V3-z ml 

where z = M 2 /m B ~ Am 2 /m 2 , and A is the gluon mass introduced to regularize the infrared divergences in the vertex 
corrections. Similar infrared divergences can also be found for the B — > Xc2 A and B — > h c K decays. The physical 
reason for the nonappearance of infrared divergences for the S-wave charmonia J ftp and n c is based on the fact that 
the virtual soft gluon in the vertex correction sees a zero total color-charge of the charm quark plus the anti-charm 
quark; whereas in the P-wave case the virtual soft gluon probes not the color charge of the charm quark and the 
anti-charm quark but the product of the color charge and the momentum of the charm quark and anti-charm quark, 
and since the total contribution from the charm quark plus anti-charm quark is nonzero (the relative momentum 
between the charm quark and anti-charm quark is essential) this leads to the appearance of infrared divergences for 
the P-wave charmonia. Experimentally, B -^XcoK has a large rate, comparable to B —tiJj'K or B -^Xc\K which are 
factorizable with the leading order V — A currents. This is certainly very puzzling in the QCD factorization approach. 

It is well known that the infrared divergences in the inclusive production and decay processes of a color-singlet 
P-wave cc state can be removed by including contributions from the higher Fock states with color-octet cc pair (say 
in S-wave) and soft gluon within the NRQCD factorization framework jE^. However, the color-octet cc pair with 
dynamical soft gluon can make contributions to the multi-body but not two-body exclusive decays. As a result, the 
infrared divergences encountered in exclusive processes involving charmonia may raise a new question to the QCD 
factorization and NRQCD factorization in B exclusive decays. Further studies are needed to seek the solution to 
remove the infrared divergences. 

As for the leading order factorizable B —>XciK decay, the calculated rate is again too small [49|. as in other leading 
order factorizable cases (e.g. B -^J/ipK and B -^n c K). 



C. B ^tp(3770)K decay 

The ^(3770) is the J PC = 1"" ID -wave charmonium with a small admixture of 2S component, and the mixing 
angle is estimated to be about 9 = —12° |5^|. Recently, Belle |54j has observed ^(3770) in the B meson decay 
B+ -> V(3770)if+ with a large branching ratio of (0.48 ± 0.11 ± 0.07) x 10~ 3 , which is comparable to B(B+ -> 
^'(3686)^+) = (0.66 ±0.06) x 10~ 3 46]. 

The B — > -0(3770)^ decay rate is calculated with QCD radiative corrections in QCD factorization [BlJ and it is 
found that the ID contribution is much smaller than 2S though both ID and 2S components in the ^(3770) wave 
function can couple to the leading order vector current and that a large mixing angle 8 = —26° is required to fit the 

observed ratio ^{^T^^, . The required large mixing angle is in contradiction with all other estimates in charmonium 

B r ( B —>ib K) 1 ° o o 

phenomenology [53|. Moreover, the calculated B — > 0(3770) K decay rate is still smaller than data by a factor of 8 
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even with the large mixing angle included. In contrast, in the B inclusive decays the D-wave charmonium states are 
expected to have large rates due to the color-octet mechanism regardless of the S-D mixing angle 55] . 

In summary, there are a number of problems in the description of B meson exclusive decays to charmonium states 
in the QCD factorization approach. The decay rates for the S-wave states such as J/ ip and r\ c are too small to accom- 
modate the experimental data. For the exclusive decays to P-wave charmonia B — ► XcjK, it is found that (except 
Xci) there are infrared divergences arising from nonfactorizable vertex corrections as well as logarithmic end-point 
singularities arising from nonfactorizable spectator interactions at leading-twist order. The infrared divergences due to 
vertex corrections will explicitly break down QCD factorization within the color-singlet model for charmonia. Unlike 
in the inclusive decays where the higher Fock states with color-octet cc pair and soft gluon can make contributions 
to remove the infrared divergences, their contributions can not be accommodated in the exclusive two body decays. 
As a result, the infrared divergences encountered in exclusive processes involving charmonia may raise a new ques- 
tion to the QCD factorization and NRQCD factorization in B exclusive decays. For the D-wave state the observed 
B —> -0(3770) AT decay rate would require a much larger S-D mixing angle than conventional expected in charmonium 
phenomenology. All these problems are likely to be related to the limitation of the color-singlet cc pair description for 
charmonium in QCD factorization. But how to solve these problems still remains an open question. New theoretical 
considerations or ingredients in QCD are needed. 

VI. PERTURB ATI VE QCD APPROACH BASED ON k T FACTORIZATION 

by CD. Lu, X.Q. Yu 

Understanding non-leptonic B meson decays is crucial for testing the standard model, and also for uncovering 
the trace of new physics. The simplest case is two-body non-leptonic B meson decays, for which Bauer, Stech and 
Wirbel (BSW) proposed the naive factorization assumption (FA) in their pioneering work (3^. The technique to 
analyze hard exclusive hadronic scattering was developed by Brodsky and Lepage |24| based on collinear factorization 
theorem in perturbative QCD (PQCD). A modified framework based on kr factorization theorem was then given in 
56. 57], and extended to exclusive B meson decays in [H |U |H . The infrared finiteness and gauge invariance 
of kx factorization theorem was shown explicitly i n | 62Tl . Using this so-called PQCD approach, we have investigated 
dynamics of non-leptonic B meson decays |34l l63l l64f . Our observations are summarized as follows: 

1. FA holds approximately for many charmless B meson decays, as our computation shows that non-factorizable 
contributions are always negligible due to the cancellation between a pair of non-factorizable diagrams. 

2. Penguin amplitudes are enhanced, as the PQCD formalism includes dynamics from the region, where the energy 
scale \i runs to yAmfe < TO fc/2, A = raj — nib being the B meson and b quark mass difference. 

3. Annihilation diagrams contribute to large short-distance strong phases through the (S + P)(S — P) penguin 
operators. 

4. The sign and magnitude of CP asymmetries in two-body non-leptonic B meson decays can be calculated, and 
we have predicted relatively large CP asymmetries in the B —* K^'ir [3J,|65| and 7T7t modes [63l 164 l6q . 

A. Formalism of PQCD Approach 

To develop the PQCD formalism for charmed B meson decays, we investigated the B — ► D^*' transition form factors 
in the large recoil region of the flW meson [6^. The B — > flW transition is more complicated than the B — > it one, 
because it involves three scales: the B meson mass mg, the flW meson mass m fl (.) , and the heavy meson and heavy 
quark mass difference, A = niB — nib ~ rn D {,) — m c of order of the QCD scale Aqcd, ™_d<») (w c ) being the D" 
meson (c quark) mass. We have postulated the hierarchy of the three scales, 

m B > m D w > A , (10) 

which allows a consistent power expansion in m D {*) /tub and in A/m D (*). 

In the B — > transition , the initial state is approximated by the bd component. The b quark decays into a c 
quark and a virtual W boson, which carries the momentum q. Since the constituents are roughly on the mass shell, 
we have the invariant masses kf ~ 0(A 2 ), i = 1 and 2, where k\ (fe) is the momentum of the spectator d quark 
m the B (£>(*)) meson. The lowest-order diagrams contributing to the B — > form factors contain a hard gluon 
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Quantities 


Cd = 0.6 


Cd = 0.8 


Cd = 1.0 


Data 


B(B° D+tt-) 


2.37 


2.74 


3.13 


3.0 ±0.4 


B(B° -» D°7T ) 


0.26 


0.25 


0.24 


0.29 ±0.05 


B(B~ -» D°tt-) 


4.96 


5.43 


5.91 


5.3 ±0.5 


Decay Modes 


C D * = 0.5 


Cd* = 0.7 


Cd* = 0.9 


Data 


B(SO ^ £>*+„■-) 


2.16 


2.51 


2.88 


2.76 ±0.21 


_» 13*%°) 


0.29 


0.28 


0.27 


0.25 ±0.07 


B(B~ -» D*°7r-) 


4.79 


5.26 


5.75 


4.60 ± 0.40 



TABLE I: Predicted B — > D^*'ty decay branching ratios in units of 10 



exchange between the b or c quark and the d quark. The d quark undergoes scattering in order to catch up with the 
c quark, forming a meson. The exchanged gluon is off-shell by 

_ fc 2 )2 „ _ _^ A 2 , (11) 

mow 

which has been identified as the characteristic scale of the hard kernels. Under Eq. 1)10(1. we have ms/m^ ^> 1, 
and the hard kernels are calculable in perturbation theory. It has been found that the applicability of PQCD to the 
B — > transition at large recoil is marginal for the physical masses tub and m D (,) |68| . 

The form factors are then expressed as the convolution of the hard kernels H with the B and meson wave 
functions in factorization theorem, 

F BDl "\q 2 )= J dx 1 dx 2 d 2 b 1 d 2 b 2 <j} B (x 1 ,b 1 )H(x 1 ,x 2 ,b 1 ,b 2 )^ DM (x 2 ,b 2 ) . (12) 

The £)(*) meson wave function contains a Sudakov factor arising from kr resummation, which sums the large double 
logarithms a s hi 2 (mBb 2 ) to all orders. The B meson wave function also contains such a Sudakov factor, whose 
effect is negligible because a B meson is dominated by soft dynamics. The hard kernels involve a Sudakov factor 
from threshold resummation, which sums the large double logarithm a s In 2 x\ or a s In 2 x 2 to all orders. This factor 
modifies the end-point behavior of the B and D^*> meson wave functions effectively, rendering them diminish faster 
in the small X\^ 2 region. 

It has been pointed out that if evaluating the hard part with a hard gluon in collinear factorization theorem, end- 
point singularities appear |69|. These singularities imply the breakdown of collinear factorization, and kx factorization 
becomes more appropriate. Once the parton transverse momenta kx are taken into account, a dynamical effect, the 
so-called Sudakov suppression, favors the configuration in which kx is not small |57j |. The end-point singularities then 
disappear. The Sudakov factor suppresses the large b region, where the quark and anti-quark are separated by a large 
transverse distance and the color shielding is not so effective. It also suppresses the x *~ 1 region, where a quark 
carries all of the meson momentum, and intends to emit real gluons in hard scattering. 



B. B — > D W 7T IN PQCD 



In this subsection we take the B — > Dir decays as an example of the PQCD analysis. In the PQCD framework 
based on kx factorization theorem, an amplitude is expressed as the convolution of hard b quark decay kernels with 
meson wave functions in both the longitudinal momentum fractions and the transverse momenta of partons. Our 
PQCD formulas are derived up to leading-order in a s , to leading power in rap/mB and in A/ mo, and to leading 
double-logarithm resummations. 

The predicted branching ratios in Table [I] are in agreement with the averaged experimental data [Til l72l l73| . 
We extract the effective parameters a\ and a 2 from our PQCD calculations. That is, our a\ and a 2 do not only 
contain the non-factorizable amplitudes as in generalized FA, but the small annihilation amplitudes, which was first 
discussed in [j4|. We obtain the ratio 1 02/01 1 ~ 0.43 with 10% uncertainty and the phase of a 2 relative to a\ about 
Arg{a 2 j a\) ~ —42°. Note that the experimental data do not fix the sign of the relative phases. The PQCD calculation 
indicates that ^4^(02/01) should be located in the fourth quadrant. It is evident that the short-distance strong phase 
from the color-suppressed non-factorizable amplitude is already sufficient to account for the isospin triangle formed 
by the B — > Dir modes. From the viewpoint of PQCD, this strong phase is of short distance, and produced from the 
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TABLE II: PQCD predictions with one angle mixing formalism (I) and two angle mixing formalism (II) and experimental data 
(in units of 1(T 4 ) of the B° -> _D ( * )0 r? (/) branching ratios. 



Decay mode 


PQCD (I) 


PQCD (II) 


Belle 


BaBar 


PDG 


B° -> D° V ' 
B° -> 5°77 
B° -> S*y 
B° -> 5*°77 


1.7 ~ 2.3 
2.4 ~ 3.0 
2.0 ~ 2.7 

2.8 ~ 3.5 


2.2 ~ 2.6 
2.6 ~ 3.2 
2.6 ~ 3.2 
3.1 ~ 3.8 


1 4+ - 6 

L - q: -0.5 

2 0+ 10 


2.41 ±0.50 


< 9.4 

< 14 



non-pinched singularity of the hard kernel. Certainly, under the current experimental and theoretical uncertainties, 
there is still room for long-distance phases from final-state interaction. 

The PQCD predictions for the B — > D*tt decay branching ratios in Table |U are also consistent with the data [ifij . 
Since mp)* and 4>p>* are slightly different from mp, and </>£>, respectively, the results are close to those for B —* Dir. 
Similarly, the ratio |a 2 /ai| and the relative phase Arg(a 2 /a 1 ) are also close to those associated with the B — > Dir 
decays. We obtain the ratio |a2/di| ~ 0.47 with 10% uncertainty and the relative phase about Argja^l^ ) ~ —41°. 
The B — > £)WjjW decays are similar to the B — > D^n° decays except the r\ -rj' mixing part The predicted 

branching ratios for these decays are shown in table [HI with uncertainties only from the rj -rj' mixing parameter 9. 

C. B s -> tv+tv- decay IN PQCD 

We also calculate the decay rate and CP asymmetry of the B s — > n + it~ decay in perturbative QCD approach 
with Sudakov resummation. Since none of the quarks in final states is the same as those of the initial B s meson, 
this decay can occur only via annihilation diagrams in the standard model. Besides the current-current operators, 
the contributions from the QCD and electroweak penguin operators are also taken into account. We find that (a) 
the branching ratio is about 4 x 10~ 7 ; (b) the penguin diagrams dominate the total contribution; and (c) the direct 
CP asymmetry is small in size: no more than 3%; but the mixing-induced CP asymmetry can be as large as ten 
percent testable in the near future LHC-b experiments and BTeV experiments at Fermilab. This small branching 
ratio, predicted in the SM, make it sensitive to possible new physics contribution. 

VII. B MESON DECAYS AND THE CKM ANGLES a, f3 AND 7 

by Z.J. Xiao 

For discussion on isospin and SU(3) relations in charmless B decays and possible new type of electroweak penguin 
or SU(3) breaking effects of strong phases, as well as large direct CP violation in charmless B decays, it is suggested 
to see ref. [2^. As a summary, we observe that in the case of SU(3) limits and also the case with SU(3) breaking 
only in amplitudes, the fitting results lead to an unexpected large ratio between two isospin amplitudes o,y 2 /ay 2 , 
which is about an order of magnitude larger than the SM prediction. The results are found to be insensitive to the 
weak phase 7. By including SU(3) breaking effects on the strong phases, one is able to obtain a consistent fit to the 
current data within the SM, which implies that the SU(3) breaking effect on strong phases may play an important 
role in understanding the observed charmless hadronic B decay modes B — > irir and nK. It is possible to test those 
breaking effects in the near future from more precise measurements of direct CP violation in B factories. 

The direct CP violation with SU(3) symmetry breaking of strong phases can be as large as follows: 

A CP (TT + n-) ~ 0.5, A C p(n°n ) ~ 0.2, A cp {tt + K~) ~ -0.1, 
A C p{-K°K ) ~ -0.1, A cp {-k°K-) ~ -0.0, A cp {-k~K°) ~ 0.1. 

For CP violation in B — + 4>Ks and possible new physics models, see ref. |2^| and references therein. 

One of the most important tasks in B factory experiments is to measure the CKM angles a, (5 and 7 the three 
inner angles of the unitarity triangle defined in the complex (p,fj) plane. At present, B — ► tttt^Kt: and many other 
interesting decays have been measured with good precision in B factory experiments. These modes together with the 
CP asymmetries of B — > J /^Ks, pp, p ±r K^ ,etc will allow us to determine the CKM angles a, (3 and 7. The isospin 
symmetry of strong interaction and SU(3) flavor symmetry are frequently used to deal with "unknown" hadronic 
matrix elements, and in some cases also plausible dynamical assumptions have to be made. 
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From a global fit done by CKM fitter group [77| the currently allowed ranges of these three angles at 95% C.L. are 

78° < a < 122°, 21° < < 27°, 38° < 7 < 80° (13) 

The angle /3 has already been well measured using the time-dependent CP asymmetries in the B® — > J/^Ks and 
related decays. The world average as given by Heavy Flavor Averaging Group [73 is 

sin 2/3 = 0.739 ±0.048, all charmonium, (14) 

= 0.692 ±0.045, all modes 

which leads to the bounds on the angle (3: 

(3 = (23.8t^)° \J (66.2+^)° , all charmonium, (15) 

= (21.9±i;?)° V (68.1+1:?)°, all modes 

The twofold ambiguity has been resolved by BaBar's new measurement of cos 2/3: a negative cos 2/3 is excluded at 
89%C.L. j2^|. We therefore believe that (3 = 24° ± 2° is a reliable measurement and can be used as an experimental 
input in the effort to determine other two angles a and 7. 

In the SM, the penguin-dominated B — > 4>K S decay should measure the same sin 2/3 as B — > J/^Ks decay. However, 
a 3.5<7 deviation from the SM prediction is observed by Belle collaboration. The most recent measurements give 



-0.96 ±0.50t^ (Belle) 
-0.47 ±0.34;°;^ (BaBar) 



wn2^ Km = i , ' ; : :»:JJ 7n\ ( 16 ) 



Although it is too early to draw any conclusion from above primary measurements, many works have been done to 
interpret such deviation as a hint of new physics contributions to the quark level b — > sss decays. 

The CKM angle a can be extracted through the CP violation in the tree-dominated B — > irir, pn and B — > pp 
decays, The B — > tt + t:~ decay currently plays most important rule in constraining a. The isospin symmetry can be 
used to eliminate the penguin pollution, but a model-independent isospin analysis of B — > 7T7t decays requires the 
knowledge of the three amplitudes A H , A 00 , A +0 and their charge conjugates A* 3 , At present, the only missing pieces 
are A 00 and A 00 . It is possible that a full isospin analysis will be done, and a extracted cleanly, by the summer of 
2005. 

Theoretically, the CP asymmetry parameter C^n and depend on a, the strong phase 5 and the ratio r = \P/T\: 

sin 2a — 2r cos S sin a 2rsin<5sina! . 

1 + r z — 2r cos o cos a 1 + r z — 2r cos o cos a 

In ref. [81|. taking 5 TT r = —0.49 ± 0.27 and (7^^ = 0.51 ± 0.19 as experimental input, we found the allowed ranges: 
(a) 76° < a < 135° for r = 0.3 ± 0.1; and (b) 117° < a < 135° and -160° < 5 < -132° if we also take the measured 
branching ratios of B — * ir + ir~, ir + ir° and K°ir + decays into account. 

In ref. by setting fx/ f-n as the SU(3) flavor symmetry breaking factor, and using the measured branching 

ratios of B — > tt + it~ and decays and the new average SV™ = —0.74 ± 0.16 and CV^ = —0.46 ± 0.13 78j as input, 
Gronau found the allowed range, a — (104± 18)° [82^. From currently available measurements of B — > decays, the 
range 19° < a < 71° is excluded at 90%C.L. From the measured time-dependent CP asymmetries in B — > /9 ± 7r =F , 
one found that the allowed range is a = 93° ± 17° (a = 102° ± 11°) if only the BaBar (Belle) measurements are taking 
into account. 

By using the SM relation a + (3 + 7 = 180°, the above bounds on a can be translated into bounds on 7 for fixed 
value of j3 = 24° ±2°. Of course, many strategies to determine angle 7 directly from B meson decays have been 
proposed recently IfE US 0> 113 > f° r example, 

• Isospin symmetry plus the B — » Kn observables allow one to constrain the angle 7. In ref. (84|, for example, 
the region 7 ~ 90° was excluded by considering a "mixed" system of B + — > ir + K° and B^ — > tt°K° decays. 



• U-spin symmetry plus the measurements of Bd — > 7T7t and B s — > ifi^, , i?d — > J/^Ks and _B S — * DjD 
B s -► Di*^^ and S d -> D^^tt^, etc. 



to determine a through the measurements of the branching ratios of six B — ► Z?/iT decays. 
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VIII. INSTANTANEOUS BETHE-SALPETER EQUATION AND ITS EXACT SOLUTIONS 



by C.H. Chang 

We propose an approach to solve a Bethe-Salpeter (BS) equation exactly without approximation if its kernel exactly 
is instantaneous, and take positronium as an example to illustrate the general features of the approach and the obtained 
solutions under it |86|. 

The key step for the approach is to start with the BS equation to derive out a set of coupled and well-determined 
integration equations for the components of the BS wave function without any approximation as a linear eigenvalue 
problem. The coupled equations may be solved analytically and/or numerically depending on the specific kernel. If 
there is no analytic solution, in principle, the coupled equations can be also solved numerically under a controlled 
(requested) accuracy. 

To solve the coupled and well-determined equations for positronium, a numerical method, which is to expand 
the equations (the eigenfunctions and the kernel) for the BS wave function components in terms of the relevant 
Schrodinger equation solutions with suitable truncation and then alternatively to solve the truncated matrix equation 
by making it diagonal, is applied, and accurate enough solutions for requests (the accuracy depends on the truncation 
mainly) are obtained. The exact solutions present precise and substantial corrections to those of the corresponding 
Schrodinger equation, which start from the order ©(w 1 ) (v is the relative velocity) for eigenfunctions, the order 0(v 2 ) 
for eigenvalues and the mixing between S (P) and D (F) components in J PC = 1~~ (J PC = 2 ++ ) states etc. 
From the lessons, we point out that such corrections are important in the effective theories such as NRQCD and 
NRQED when we consider the relativistic corrections 0(f) ~ O(jj). Namely one cannot 'forget' the same order 
corrections involved in the solutions of the bound states accordingly. 

Moreover, we also point out that there is a questionable step in the classical derivations for an instantaneous BS 
equation, if one is pursuing the exact solutions without approximation. Namely it has confused the differences between 
the instantaneous BS equation and Breit equation (more precise illustration in |87|l. 



IX. NEW PHYSICS EFFECTS IN B DECAYS 

by X.G. He and C.S. Huang and Z.J. Xiao 



A. Time Dependent CP Violation in B -» K*j in SUSY 

In this work we study time dependent CP asymmetries in B — > K*"f — > ttKs"/ and B — > <f)Ks using constraints from 
B — > Xsj in SUSY model with large Left-Right (or Right-Left) squark mixing induced gluonic dipole interaction. 
Since this work was reported at the meeting, new results on B decays have been reported |88, 89, 9(J,|9l|], in particular 
the direct CP asymmetry Acp in B° — + K~ir + of 88] — 0.114 ± 0.020 measured by Babar and Belle, we include the 
constraint from Acp{K~tv + ) in our study [92|. 

In the SM, the Hamiltonian for the B decays to be considered is well known which is of the form [9j| , 

„ 12 

H = -^[V ub V: s { Cl O x + c 2 2 ) - J2 VjbVfAOi], (18) 

where Vy are the CKM matrix elements. Cj are the Wilson coefficients for the operators Oi which have been evaluated 
in different schemes. Values from NDR scheme will be used |93|. We will not display the full sets of Oi and Ci here, 
but only give the definitions of the gluonic and photonic dipole operators On and 0\i for the convenience of later 
discussions. They are given by 

On = /^sa^Grr Q K(l + 75) + m s (l-7 5 )]6, 

°" = ^sa^F^[m b (l+75)+m 5 .(l- 75 )]6, (19) 

where T a is the color SU( 3) g enerator normalized to Tr(T a T b ) — S ab /2. G M „ and F^ v are the gluon and photon field 
strengths. In the SM [HIMlil dx = -0.151 and c 12 = -0.318. 

When going beyond the SM, there are modifications. In SUSY models, exchanges ofgluino and squark with Left- 
Right squark mixing can generate a large contribution to at one loop level 961 197| since their interactions are 
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strong couplings in strength and also enhanced by a factor of the ratio of gluino mass to the b quark mass |98| ■ We 
will concentrate on the effects of this interaction. In general exchange of squarks and gluinos can generate non-zero 
Cn,i2 for dipole operators with 1 + 75, as well as with non-zero c' n 12 for dipole operators with 1 — 75. 

The Wilson coefficient c^ 8 ^ for SUSY contribution obtained in the mass insertion approximation is given by, for 
the case with 1 + 75 [13 , 



=11 i m g) 



%/27ra s (mg) 6 b L s R m g 
G F m\ VtbVts m b 



>n<: 



12 



Lrpm- VtbVts m b 



G (Xgq), 



•g • <.u'ts 

-2 1 (Mi^ ™2 



G (.t) = 
F (x) = - 



;[22 - 20x - 2a; 2 + (16a; - a; 2 + 9) ln(x)] 
3(1 -a;) 4 : 
4a; [1 + 4a; - 5a; 2 + (4a; + 2a; 2 ) ln(x)] 
9(1 -a;) 4 ' 



(20) 



where 8f^ R parameterizes the mixing of left and right squarks, x gq = rn^/m- is the ratio of gluino mass m g and squark 
mass rriq. The Wilson coefficients c/^f for the case with 1 — 75 can be easily obtained by replacing the Left-Right 
mixing parameter S b p R by the Right-Left mixing parameter S RL . 

From the expressions in ea.l|2()[l. one can see that the SUSY contributions are proportional to 



If m-g is of order 

a few hundred GeV, there is an enhancement factor of (mg/mb)(m^y /rn~) for the SUSY dipole interactions. In this 
case even a small 8°f R RL , which can easily satisfy constraints from B — 
effect on rare B decays. 

We first consider constraint on the SUSY parameters 8 b p R RL from B 

has been measured to a good precision with [8£j (3.54+Q2g) x 1CU 4 . Although experimentally CP asymmetry in 
B — » Y57 has not been established, there are constraints from experiments with 0.005 ± 0.036. We follow ref.|95l| 
to carry out the evaluation of the branching ratio and CP asymmetry and use them to constrain the parameters. 

i g = 300GeV and rriq in 



B mixing and other data, can have a large 
-> Xsj- The branching ratio of this process 



In Figure [5] we show the allowed ranges for the absolute value of 8^ R and its phase r for 



the range 100 ~ lOOOGeV at the one a level. We find that the constraints from BR(B — > Y57) and Acp{B — > X$"f) 
are similar. Using the allowed parameters, one can obtain the allowed en through eg 1201 and to study implications 
for other rare B decays. 




j 1 1 1 j 1 i_ 
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FIG. 2: The one a allowed ranges for the SUSY parameters \8 h j^ RRL \ and the phase r taking mg = 300GeV and rriq in the 
range 100 ~ 1000 GeV. The light-dark doted areas are the allowed parameter spaces from BR(B — > X 3 ^y) and Acp(B — > Xs"f) 
constraints. The dark doted areas are allowed ranges by Acp(B° — > K~tt + ) constraint. Figure la (on the left) and Figure lb 
(on the right) are for the dipole operators with 1 + 75 and 1 — 75, respectively. 

The recently measured CP asymmetry Acp(B° — > K~ir + ) can be reproduced and the SUSY parameters can be 
further constrained. We follow ref.[9^| to calculate the SUSY dipole contribution to B — * Ktt. In our numerical 
analysis we take the CKM parameters to be known, with the standard parametrization s±2 = 0.2243, S23 = 0.0413, 
S13 = 0.0037, 813 — 1.05, which is the central value given by the Particle Data Group ^(|- With the SM amplitudes 
obtained and the default values for the hadronic parameters used in Ref. [9^, we obtain the CP asymmetry Acp{B° — > 
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K~ir + ) in the SM to be 0.15. This is different in sign with the experimental value. When SUSY dipole interactions are 
included the experimental value can be reproduced. For example with nig = — 300GeV, S^p — 2.62 x 10~ 3 e 23Sl , 
5 RL = 4.31 x io- 3 e 1 007i the asymmetry A C p(B° — > K~tt + ) is approximately -0.114. With the same set of SUSY 
parameters, we have Br(B — + Xgj) — 3.48 x 10 -4 , Acp(B — > Xsj) = 0.016. It is clear that the CP asymmetry 
Acp(B° — > K~tt + ) can be brought to be in agreement with data at one sigma level when SUSY gluonic dipole 
interactions are included. 

To see how the CP asymmetry provides stringent constraint on the SUSY flavor changing parameters, we show in 
Figure|21the parameter space allowed from Acp(B° — > K~n + ) (the dark doted areas) on top of the allowed ranges by 
B — > Y57 constraint alone at the one a level. We see that the CP asymmetry in B° — > K~tt + considerably reduces 
the allowed parameter space. 

We now study time dependent CP asymmetries in B — ► K*^ — > TT°Kg^ and B — ► 0if,g. There are two CP violating 
parameters Af and Sf which can be measured in time dependent decays of B and B produced at e + e~ colliders at 
the T(45 l ) resonance, A CP (t) = A s cos(AtAm B ) + S } sm(AtAm B ). 

For B° -> if*°7 -> 7r ifs 7 and B° -> K*°j -> 7r Ks7, we obtain pi Iflsl fTcTnl| 



c _ r, Im [(QB/pB){ci2c' 12 )} 

bK "< | C 12| 2 + |ci 2 | 2 • ( ' 

To the leading order Apc* 7 is the same as Acp(B — > X57). Note that the hadronic matrix element < X* |scr A " y (l ± 
7s)6|-B > does not appear making the calculation simple and reliable. In order to have a non-zero Sr-* 7 both C12 and 
c' 12 cannot be zero. 

In the SM the asymmetries A K . 1 and S K * 7 are predicted to be small with [il EH3 ^|^ 7 » 0.5%, « 3%. 

With SUSY gluonic dipole interaction, the predictions for these CP asymmetries can be changed dramatically [98). 
With the constraints obtained previously, we find that the parameter qs/pB is not affected very much compared with 
the SM calculation. To a good approximation qb/pb = e~ 2i/3 . 

A Large gluonic dipole interaction also has a big impact on B — > (f>Ks decays |101[ . In the SM, A$k s is predicted 
to be very small and S$k s is predicted to be the same as Sj/^k s = sin(2/3). With SUSY gluonic dipole contribution, 
the decay amplitude for B — » 4)Kq will be changed and the predicted value for both A^Ks anc i S<j>K s can be very 
different from those in the SM |101| . We again use QCD factorization |9?| to evaluate the amplitude. 




FIG. 3: The allowed time dependent CP asymmetries in B° K*j — > Ksn°j and B° -v 4>Ks. 

The results are shown in Figure 3. The current values of Sr-* 7 and A^« 7 from Babar (Belle) are |102| : 0.57 ± 0.32 ± 
009(-0.00 ± 0.38), 0.25 ± 0.63 ± 0.14(— 0.79±^Jg ± 0.09), respectively. From Figure 3, we see that the allowed ranges 
can cover the central values of Sk*~/ from Babar and Bell, but it is not possible to obtain the central value of Ak* 7 
by B elle. Future improved data can further restrict the parameter space. Both Babar and Belle have also measured 
|j04| A C p(B- -> K*~i) with ranges -0.074 - 0.049(Babar) and -0.015 ± 0.044 ± 0.012(Belle). In the model we 
are considering, the CP asymmetries Ak-* 7 and Acp(B~ — > K*~j) are the same. The results for the charged B CP 
asymmetry is consistent with data. 

The time dependent asymmetry in B — > 4>Ks is a very good test of CP viol ation in the SM. Experimental measure- 
ments have not converged with the current values of Babar (Belle) given by [U|l03| 0.00±0.23±0.05(0.08±0.22±0.09), 
and 0.50±0.251o;q4(0.06±0.33±0.09) f° r ^<t>K s an d SipK s > respectively. These values are considerably different than 
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the value reported by Belle last year of [Toll S+Ks = -0.96±0.50±g;??. From Fi gure 3 we see that the current data of 
A^Ks an d S^Ks can be easily accommodated by the allowed ranges. We also note that the allowed ranges can cover 
last year's Belle data. Since the error bars on the data are large, no definitive conclusions can be drawn at present. 

B. CP asymmetries in B -> (j>Ks and B -> rj'Ks in MSSM 

The time dependent CP asymmetry in B — ► <j>Ks 

S<j,Ks = -0.39 ± 0.41, 2002 World - average 

S<j>K s = -0.15 ± 0.33, 2003 World - average (22) 

is especially interesting since it deviates greatly from the SM expectation 

S^k s = sm{2f3(q>K s )) = sm(2f3(J/^K s )) + 0(\*) (23) 

where A ~ 0.2 appears in Wolfenstein's parametrization of the CKM matrix. If the error in the average is inflated 
by the scale factor, the 2003 World-average will be — 0.15± 0.69, which deviates from the SM only by 1.3c, i.e., the 
statistical basis of the effect is weak at present. Although the impact of these experimental results on the validity of 
CKM and SM is currently statistically limited, they have attracted much interest in searching for new physics and 
it has been shown that the deviation can be understood without contr adict i ng t he smallness of the SUSY effect on 
B — ► J/ipKs in the minimal supersymmetric standard model (MSSM) |l06l Il07| . 

Another experimental result is of the time dependent CP asymmetry S^Ks m B ~ * v' Ks 

S^Ks = 0.02 ±0.34 ±0.03 BaBar 

= 0.43 ±0.27 ±0.05 Belle (24) 

which deviates sizably from the SM expectation. Although both the asymmetries S$k an d S^k are smaller than the 
SM value, Sqk is negative and S^k is positive, which implies that new CP violating physics affects B — > 4>Ks in a 
dramatic way but gives B — > rj'Ks a relatively small effect. Because the quark subprocess b — > sss contributes to both 
B — > <pKs and B — > rj'Ks decays one sho uld simultaneously explain the experimental data in a model with the same 
para mete rs. It has been done in Ref. |108| in a model- independent way in the supersymmetric (SUSY) framework. In 
Ref. 108] the analysis is carried out using the naive factorization to calculate hadronic matrix ele men ts of operators 
and the neutral Higgs boson (NHB) contributions are not included. As we have shown in a letter |l07j that both the 
branching ratio (Br) and CP asymmetry are significantly dependent of the a s corrections of hadronic m atrix elements 
and NHB contributions are important in MSSM with middle and large tan/3 (say, > 8). In this paper 109] we study 
the B — > 4>Ks and B — ► rj'Ks decays in MSSM by calculating hadronic matrix elements of operators with QCD 
factorization approach and including neutral Higgs boson (NHB) contributions. We calculate the Wilson coefficients 
of operators including the new operators which are induced by NHB penguins at LO using the MIA with double 
insertions. We calculate the a s order hadronic matrix elements of the new operators for B — * <pK s and B — ► rj'K s . 
We analyze constraints from relevant experimental data among which the new CDF bound of B s — > is used. It 

is shown that the recent experimental results on the time-dependent CP asymmetries in B — > <j)Ks and B — > rf ' K$ 1 
S^k is negative and S^k is positive but smaller than 0.7, which can not be explained in SM, can be explained in 
MSSM if there are flavor non-diagonal squark mass matrix elements of 2nd and 3rd generations whose size satisfies 
all relevant constraints from known experiments (B — > Xs"/,B S fi + [i~,B — > A s /i + fi~ , B — -> l s j, AM S , etc.). In 
particular, we find that one can explain the experimental results with a flavor non-diagonal mass insertion of chirality 
LL or LR when a s corrections of hadronic matrix elements of operators are included, in contrast with the claim in 
the literature. At the same time, the branching ratios for the two decays can also be in agreement with experimental 
measurements. 



C. B decays in model III and TC models 

In Ref.[ll0| the authors made a systematic study for the new physics contributions to the charmless two-body 
hadronic decays of B and B s meson induced by the neutral- and charged-Higgs penguin diagrams in model III: the 
third type of the two-Higgs-doublet models. They found that the new physics enhancements to the branching ratios 
of those penguin dominated B meson decay modes can be significant, about 30% to 50%, but the corresponding new 
physics corrections to CP violating asymmetries are generally small. Also in model III, the penguin contributions to 



18 



the de cays of b — > sg, B — > X s -f , X s ^fj, V-f (V — K*,p, uj), l + l and B° — B° mixing have been studied, for example, 
in Ref.[m| 

In the framework of various technicolor models, many B meson decay modes, such as the radiative decays B —* X s "f 
and B ~ + A s 77, the semileptonic decays B — > K^*H + l~ and B — > A s i/z/, the charmless two-body hadronic decays 
B — ► /ii/i.2) have been calculated recently |l!2| . Strong constraints on the parameter space of technicolor models have 
been found from these studies. 



X. SUMMARY 



In this B physics workshop, a number of interesting topic of B physics is discussed. The heavy quark effective theory 
and heavy quark effective field theory are intensively discussed and compared. The two approaches for hadronic B 
decays, PQCD and QCD factorization are also discussed and compared. A lot of questions are raised. The hadronic 
wave function including its twist 3 part are derived. The contribution of new physics especially super symmetry 
are discussed. Most scientists and students benefit a lot from the lectures and discussions. For a summary, Prof. 
Yue-Liang Wu writes In his transparencies: 

• BOTTOM IS BEAUTY ! b <-> b 

• BEAUTY FROM TOP ! b <- t 

• BEAUTY TO UP ! b -> u 

• BEAUTY TO CHARM ! b -> c 

• BEAUTY TO STRANGE ! b -> s 

• PHYSICS REMAINS AT BOTTOM , LCQCD /HQEFT /QCDF /PQCD IS VERY USEFUL 

• MORE NEW PHYSICS BEYOND SM, PHYSICS IS BEAUTY & CHARM 
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